function [ lijst ] = relationalProfile( house,sens1,sens2 )
    
    start_time=house.ss(:,1);
    sensor_id=house.ss(:,3);
    data=[start_time sensor_id];

    sensorfilter=(data(:,2)==sens1);
    data1=data.*repmat(sensorfilter,1,2);   data1=data1(any(data1,2),any(data1,1));
    sensorfilter=(data(:,2)==sens2);
    data2=data.*repmat(sensorfilter,1,2);   data2=data2(any(data2,2),any(data2,1));

    lijst=[];
    for i=1:length(data1)
        shortest=1;
        for j=1:length(data2)
            inbetween=data1(i,1)-data2(j,1);
            if(abs(inbetween)<1 & abs(inbetween)<abs(shortest))
                shortest=inbetween;
            end
        end
        lijst=[lijst;shortest];
    end

    %scatter(lijst,lijst./lijst);
    %X=-1.5:0.000001:1.5;    plot(X,normpdf(X,mean(lijst),var(lijst)));
    
end

